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The 'present paper critically discusses the method of detection^ the magnit'ude and the 
rate of occurrence of sudden disturbances in the motions of some short-period comets. 
The disturbances have recently been suggested as potential indicators of collisions be- 
tween the comets and interplanetary boulders — minor objects whose existence was pre- 
dicted by M, Harwit in 1967. The character of explosive phenomena, caused by an 
impact of such a boulder on a comet's nuclear surface, depends significantly on the 
surface texture of the target body. To advance our uj^rstanding of the impact mechcv- 
nism., a method is suggested which wouM_supply a good deal of the missiy infomuUton 
about the structure and optical properties of nuclear surfaces from precise photometric 
observations of cometary nuclei at large solar distances. 


THE IMPACT HYPOTHESIS 

R ecent extensive dynamtcali studies of a 
number of short-period comets by Marsden 
(1969, 1970), by Yeomans (1972), and by 
Marsden and Sekanina (1971) resulted in a 
discovery of easily detectable disturbances in the 
motions of the comets we call ^erratic^ P/Biela, 

P /Brorsen, P /Giacobini-Zinner, P /Perrine-Mrkos, 
P/Schaumasse, possibly also P/Forbes and 
P/Honda-Mrkos-Pajdu§akov4u The disturbances 
differ from the regular nongravitational effects 
and seem to take form of sudden impulses of 
about 1 m/s, perhaps preferably at larger solar 
distances. Their interpretations in terms of 
processes stimulated by internal cometary sources 
of energy have been discarded on various grounds. 

In contrast, hjiiervelocity impacts of small objects 
seem to be consistent with empirical evidence. 

The existence of interplanetary boulders has been 
predicted by Harwit (1967), of spatial density as 
high as 10“*^ g/cml Application of the mechanism 
of crater formation at hypervelocity impacts 
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suggests that the observed disturbances can be 
generated by collisions of the boulders with low 
density comet nuclei, if the comet-to-boulder 
mass ratio is about 10®. As a result of such an 
impact the comet would lose as much as 10 
percent of its mass. Repeated impacts can easily 
result in a splitting of the nucleus, or its complete 
disintegration in a relatively short period of time. 
With HarwiPs space mass density of boulders the 
proposed hypothesis predicts an average rate of 
some five impacts per 100 revolutions for a comet 
1 km in diameter. To produce an impulse of 1 m/s 
the average boulder should be 10® g in mass, or 
3 to 10 m in diameter, depending on its mass 
density; the comet would be 10^^ g in mass, and 
0.2 g/cm^ in density. Such a nucleus can be com- 
posed of snows mixed with highly porous dust 
grains. Impacts of the same boulders would not 
measurably affect motions of the comets with 
heavy compact cores like P/Encke. 

These have been some of the main conclusions 
formulated by Marsden and Sekanina (1971) 
from their extensive study of the motions of the 
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^erratic^ comets. In the present paper we discuss 
numerical values of some of the fundamental 
constants of the problem in greater detail. We 
will refer to the above paper as to Paper 1, 


RADIAL AND TRANSVERSE COMPONENTS 
OF THE NONGRAVITATIONAL EFFECTS 
IN THE DAILY MEAN MOTION 

An acceleration, the radial component of which 
is Zi (positive outward from the Sun) and trans- 
verse component Z 2 (perpendicular to Zi in the 
orbit plane, positive in the direction of motion), 
applied at a solar distance r generates an in- 
stantaneous rate of change in the daily mean 
motion 




3c sin y „ 3(1 — 

a(l-e*)‘« r ’ 


( 1 ) 


where v is the true anomaly, a and e are the semi- 
major axis and eccentricity of the orbit. We 
accept that the nongravitational acceleration 
varies with the solar distance r, 

Zi = Aif{r) 


hand, from equation (2), 

A = -Axkhp-^i^ (6) 

SO that 

( A^)r = - {2ir)~^AJcha-^i^ (7) 

It is convenient to convert Afx of equation (4) 
to AT 2 and (A>i)r of equation (7) to ATi, the 
effective rates of delay (AT,- >0) or advance 
(AT, <0) in the perihelion passage per revolution 
due to, respectively, the transverse and radial 
components of the acceleration of equation (2) . 
Expressing AT^ in days per revolution per revolu- 
tion we find 

ATi = 58Ai/2a3/2 

ATa= 

The integrals Ii and 1 2 are of the same order of 
magnitude. For a typical short-period comet ratio 
AT 2 /AT 1 is about 10, even when the radial com- 
ponent is almost an order of magnitude larger 
than the transverse component. 

DYNAMICAL DISTURBANCES 


/(1AU) = 1 (2) 

x = l,2,3 

where A i arc the acceleration components at 1 AU 
in units of the solar gravitational acc(^leration at 
1 AU. Upon integrating over a revolution period, 
and writing 


Is= 



(3) 


we obtain the change in the daily mean motion 
per revolution : 

Ajjl= —3A2kIia~^^^ (4) 

There is no contribution from the periodic varia- 
tions in the radial component of the acceleration. 
However, there is a secular effect from Zi because 
of its contribution to the “effective” gravitational 
constant. On the one hand we have 



= 27TjP Ak — 2irp'^^^^a^^^(A}i)r (5) 


where p = a{l—e^) and (Aju)r is the change, per 
revolution, in the dail}^ mean motion due to the 
change in the Gaussian constant k. On the other 


The computer programs used by the authors, 
mentioned earlier for calculating the orbital 
elements and nongravitational parameters from 
comet observations, are designed to search for 
smoothly, continuously varying deviations from 
the gravitational law. If a disturbance is detected 
by the program in a comet's motion, contradicting 
the above assumption, the integration procedure 
does not necessarily fail. What does happen 
depends much on the number of the comet's 
apparitions linked. A solution may be found, 
which gives quite an acceptable distribution 
of residuals, but the nongravitational parameters 
are inconsistent with those computed from the 
comet's adjacent apparitions not including the 
dynamical anomaly. Figure 1 shows an example 
of such a forced solution. The regular nongravita- 
tional effects shape the continuous background 
SABEF Z. Bctwxen Ib and is a disturbance BCTDE 
is superposed on the quiescent phase BE. If an 
attempt is made to link apparitions between Ia 
and tzj and the nongravitational effects are 
allow'ed to vary exponentially with time, 

A^exp{ — Bt) ( 9 ) 

the empirical fit yields the curve A'CDFZ' such 
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UioTTBE 1. — Sudden dynamical anomaly, or disturbance, 
interferes with continuous, quiescent-phase nongra\i- 
tational effects in a comet’s motion. The disturbance 
can be detected by means of a ‘forced’ solution: the 
cross-shaded areas compensate the one-way shaded 
areas. 

that 

area(AA'CBA) + arca(DFED) 

= area(CTDC) + area(FZZT) (10) 

The coefficient B of the secular variations — 
negative in the quiescent phase in figure 1 
suddenly becomes positive due to the disturbance* 
Thus, figure 1 is a very obvious demonstration 
that whenever a disturbance is involved, B comes 
out fictitious. It is easy to understand that the 
sign and magnitude of B depends not only on the 
disturbance-to-background ratio, but also on the 
selected span of time. For example, B would be 
strongly’' negative, if we tried to link apparitions 
between ts and Ie- Moreover, the forces generating 
the quicscent-phasc and disturbance effects may 
work in opposite directions, and we may fail to 
find a satisfactory solution of the form of equation 
(9) and must accept another empirical form. If 
only three apparitions are linked it is always 
possible to find a satisfactory solution with 
constant A. Figure 1 corresponding to this case 
would have a staircase shape, and the general 
rule, equation (10), would again be in power. In 
practice, however, the validity of equation (10) 
is only approximate. The reason comes from the 
difference between the real orbit (with the un- 
known profile of the disturbance) and the fictitious 
orbit found by the forced continuous solution. The 
differential perturbations, predominantly due to 


Jupiter, along the two orbits should be taken into 
account in equation (10), In practice, the per- 
turbations are very^ small unless the comet makes 
a close approach to Jupiter during the critical 
period of time. Unfortunately, these encounters 
are fairly frequent and often limit our results in 
accuracy. 

SUDDEN IMPULSES 

We do not— and practically cannot— have direct 
evidence of the character of the dymamical dis- 
turbances affecting the "erratic' comets. We guess 
that they take form of discrete discontinuities 
(sec Paper 1), because so far it has alvay^s been 
found that observations from only- the minimum 
number of apparitions, necessary for the least- 
squares procedure to work, can satisfactorily’^ be 
fitted whenever a disturbance is involved. Outside 
that span the forced solution completely fails. 
Typically, there are long intervals of quiescent 
phase before such a disturbance (P/Giacobini- 
Zinner) , or after it (P/Biela) , or both before and 
after (P/Schaumasse) . 

If wc take the disturbance in the form of a 
sudden impulse iu fig- 1) 

estimate the quiescent-phase background, we can 
determine the impulsive increment in the orbital 
velocity, AUr, associated \rith the discontinuity^ 
from the difference between the disturbed and 
quiescent nongravitational parameters. The im- 
pulse corresponds to the area BCTDEB in 
figure 1. 

The component of the nongravitational accelera- 
tion along the orbital velocity vector, Zp, is 
given by 

Z,= (ZiFi+Z2F2)U-^ (11) 

where Zi, Z 2 are identical with those of equation 
(2), V is the orbital velocity, 7i and F 2 its radial 
and transverse components respectively. In- 
tegrating over the revolution period, we have 
from equation (11) 

fp (2 

^p= j Z, (U=kA2 J f(r)r dv (12) 

where v is the true anomaly; k = 29.8X1(F, if yp is 
to be given in meters per second. The expression, 
equation (12) , is independent of Zx for the reasons 
discussed above. 

Let lAquiesc and ^autrb be lA for the quiescent phase 
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and the disturbed period, respectively. The 
impulsive increment in the orbital velocity 
associated with the disturbance is then 

V (l^distrb ~“^quiesc) ^ 

+ (effect of differential perturbations) (13) 

where v is the number of revolutions covered by 
the forced fit. The positive AT» means the comet is 
effectively decelerated, the negative means accel- 
erated. If A 2 has been allowed to be subject to 
secular variations, its effective value during the 
period of time covered by the forced fit must be 
used in equation (12). 

The total impulsive velocity, AT, associated 
with the disturbance cannot be derived from its 
component along the orbital velocity vector, 
because the angle between the impulse and the 
orbit tangent is not known. Assuming that the 
discontinuity in motion is due to a collision with 
a small object moving in a circular orbit around 
the Sun in the comet’s orbital plane, we have 
derived in Paper 1 the following formula for the 
mean quadratic relative velocity between the two 
colliding bodies, averaged by integration over a 
revolution period : 

w= (iy2)i/2 

[ 4 p/2 

(14) 

where K.{7n} is the complete elliptic integral of 
the first kind with modulus m. Similarly we can 
calculate the mean quadratic component of the 
relative collision velocity in the direction of the 


Table 1 . — Ratio aT/AV» vs Eccentricity^ c 


€ 

AF/A7, 

0.1 

2.24 

0.2 

2.22 

0.3 

2.19 

0.4 

2.15 

0.5 

2.10 

0.6 

2.04 

0.7 

1.96 

0.8 

1.87 

0.9 

1,75 


comet’s motion: 

L 7T 


" 11/2 




( 15 ) 


Relating the direction of the impulsive velocity 
to that of the relative collision velocity w^e now^ 
can find AT averaged over a revolution period: 

AT=AF. 


X 


1 - 


l-(l-e2) 


1/2 


1 - 1/2 


2 - - {l-ey^^K[{2ey^^{l+e)~ 




( 16 ) 

The ratio AF/AT,, is listed in table 1 as a func- 
tion of the eccentricity. 


IMPULSIVE VELOCITIES ASSOCIATED 
WITH THE DISTURBANCES OBSERVED IN 
MOTIONS OF THE 'ERRATIC’ COMETS 

The method described in the preceding section 
can give a reasonable estimate of the impulsive 
velocity, particularly if there is no close approach 
to Jupiter involved. Table 2 lists AF, and AF 
obtained in this way for the 'erratic’ comets and 
compares them with the values derived by 
Marsden for P/Schaumasse, P/Perrine-^Irkos 
and P/Biela, and by Yeomans for P/Giacobini- 
Zinner, w^ho have used a different approach. 
These authors have computed what wq call AF^; 
from the difference betw^een the observed time of 
perihelion passage and that extrapolated from a 
quiescent phase. 

The fundamental difference betw^een the tw'o 
methods is that the one w^e suggest tends to 
smooth the disturbance out and represents there- 
fore a low^er limit of the most probable impulsive 
velocity. On the other hand, the method applied 
by Marsden and by Yeomans extrapolates, and 
therefore tends to exaggerate the effect of the 
disturbance. Indeed, table 2 clearly shows that 
our values of the impulsive velocity are system- 
atically smaller. In any case, the table suggests 
that 1 m/s, w'hich \vas accepted in Paper 1 for the 
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Table 2. — ‘Erratic' Cornels: Impulsive Velocities Associaied With Dynamical Disturbances 


Comet “ 

Disturbance 
between — 

P/Biela 

1772/1805 
<* 1842/43 

P /Brorsen 

' 1846/73 
1873/79 

P/Schaumassc 

1927/43 

P /Perrine-M rkos 

1955/68 

P/Giacobini-Zinner 

1959/65 


aF„ 

from eq. 
(13) (m/s) 

aF 

from eq. 
(16) (m^) 

AF. 

derived 

otherwise^ 

(m/s) 

Comet’s mean motion 
effectively® 

1 

In quiescent 
phase 

By 

disturbance 

0.8 

1.5 


A 

D 

0.1 

0.2 

« 1 3/2 

A 

D 

(0.3) 

(0.5) 


A 

7 

1.6 

3.0 


A 

D 

0.2 

0-5 

1.9 M, 

A 

A 

1.5 

3.0 

3 . 5 M 1 

A 

A 

0.6 

1,1 

1.4 Y 

D 

D 


• Data on P/Honda-Mrkos-Pajdu-dkovd and P/Forbcs arc inconclusive. 

>> 3/i = Marsden (1970); Afi = Marsden (1971) [also in Marsden, Sekanina (1971)]; F = Veomans (1971). 
® A = accelerated; D = decelerated. 

^ Satellite nucleus at splitting. 

® Velocity of separation essentially in sunward direction, 
f Impulsive velocity rather uncertain. 


typical impulsive velocity, seems indeed to be a 
representative value for the * erratic^ comets. 

RATE OF DISTURBANCES 

The rate of dynamical disturbances is another 
critical quantity for the impact hypothesis. 
Table 3 lists the number of observed disturbances, 
0fl, for the 'erratic’ comets; the length, of their 
trajectories swept out between the first and last 
observed apparitions of each of the comets; and 
the disturbance rate, r, defined as the number of 
disturbances per 100 AU : 

r=100Dl/^ (17) 

For the kno\\m 'erratic’ comets the average 
observed rate is about one disturbance per ITO AU 
(see column 6 of table 3), a value about four 
times as high as the one predicted from the impact 
hypothesis in Paper 1 . However, it is easy to show 
that the data of table 3 are strongly affected by 
observational selection. 

The number of boulder impacts on a cometary 
nucleus is statistically proportional to the volume 
of space swept out by" the comet. The volume is 
given as a product of the comet’s collisional cross- 
section and the length of its trajectory. Ideally, 


91 should depend linearly on t in figure 2. However, 
since we deal with observed lengths of trajectories 
that are very" short compared to the rate of pre- 
sumed impacts, statistical dispersion is significant 
and the observed disturbance rates for individual 
comets differ widely^ from each other. We identify 
only the best observed 'erratic’ comets. It is 
therefore logical that the upper left comer of 
figure 2, the area of the highest disturbance rates, 
is populated most, whereas the strongest bias 
takes place along 91 = 0. Observational selection is 
also responsible for a factor of three beUveen the 
mean least-square rate 91// (dot-and-dashed line) 
and the differential rate d^/dt (dashed line) . It is 
the latter that should more properly match the 
unbiased disturbance rate. Indeed, the rate pre- 
dicted from the impact hypothesis (solid line) 
agrees with the slope d9l/ dt. 

For the above reasons it is convenient to write 
the number of the 'erratic’ comets with dis- 
turbance rates between P and r-(-dr in the form 

dNc{V)^NM^)dvJ 0(F) dP = l (18) 
*'0 

and to see how the total number, of the 'erratic’ 
comets among known short-period comets and 
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Table 3. — ‘Erratic^ Comets: Observed Rates of Dynamical Disturbances 


Comet 

First/last 

apparition 

Number of 
revolutions 
covered 

Length of 
swept-out 
trajectory 

(AU) 

Number of 
detected 
disturbances 

Disturbance 

rate 

(per 100 AU) 

P/Biela 

1772/1852 

12 

226 

2 

0.88 

P/Brorsen 

1846/79 

6 

95 

2 

2.10 

P/Schaiimassc 

1911/60 

6 

131 

1 

0.76 

P/Perrine-Mrkos 

1896/1968 

11 

210 

1 

0.48 

P/Giacobini-Zinnor 

1900/66 

10 

187 

1 

0.53 

P/Honda-Mrkos-Pajdiis4kov^ 

1948/69 

4 

61 

1? 

M.23 

P/Forbes 

1929/61 

5 

100 

1? 

•0.75 


Tho disturbance rate of this comet is weighted by a factor of ^<4 to allow for the uncertainty as to whether the dis- 
turbance indeed occurred. 



Figure 2. — Obsen^cd dynamical disturbances. Solid cir- 
cles: definite 'erratic' comets. Open circles: probable 
'erratic' comets. Solid line: disturbance rate predicted 
from the impact hypothesis (see Paper 1). Dot-and- 
dashed line: least-square solution to the mean ob- 
serv'cd disturbance rate, Dashed line: least- 

square solution to the mean-differential disturbance 
rate, 

the characteristic disturbance rate, given by 

p 

<^(r) rfr= f .^>(r) rfr (19) 

•'Fc 

depend oil the choice of the distribution function 
0(r). The number of comets with disturbance 
rates higher than F, 

Nc^{T)=Nc f dy (20) 

•'r 

proves the most useful quantity for practical 
trials, because our statistics of disturbances is 
relatively complete for very high values of F. 

Approximating 0(F) first by the ]\Iaxwcllian 


velocity distribution function, with T„, being the 
most frequent F, 

0(F) dF = 4r-^^2p2p^-3 exp[~ (F/F,n)2] dV (21) 
we find 

A.+ (F)=A,{2ir-^/2rr^-i 

X exp[- ( F/F,n) 2]+erfc ( F/F^) } (22) 

where 

erfc(x) = 1 — 27 t “*/2 f exp( — a:^) dx (23) 
•'o 

The characteristic rate is 

T,= h0S7r^ (24) 

The fit to the empirical data is, however, un- 
satisfactory. 

Alternatively, we can assume that 0(F) has 
the form of a two-dimensional ]\Iaxwellian velocity 
distribution. This assumption seems to be more 
plausible in view of presumably low obliquities 
between orbital planes of the short-period comets 
and interplanetary boulders. Then 

0(F) dT^TTm-^ exp[- r/2TJ'] dV (25) 

A,+ (F) =Nc exp[-FV2F,„2] (26) 

and 

Fe=F,n(l0ge 4)1/2 (27) 

There is an improvement upon equation (22) in 
matching the data of table 3, but the fit is still 
poor for Ac"^<3. 
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Let us next accept an exponentially decreasing 
distribution 

4>{T) exp[-r//3] dr (28) 

so that 

Nc^ {T)^Nc exp[~ r/^] (29) 

and 

r. = ^-loge2 (30) 

The fit is now good except for Nc^ =1. 

Finally, if we take 

(^(D dr=(2r)“ir-i/2exp[-r^^Vf]dr (3i) 

at +( r) =Nc exp[-- (32) 

r.= (Mog.2)2 (33) 

we get a very good fit to all the seven data. 

Unfortunately, the testing data are too scanty 
to resolve the ambiguity in <f>{T) unequivocally, 
and the best fit does not necessarily mean the best 
solution. Indeed, a fit at least as good as that by 
equation (32) is obtained from 

Ar + (r)=2.7r-i'* (34) 

which gives no prediction for Nc whatsoever. 

Table 4 lists the characteristic rate Tc and the 
extrapolated Nc for the four applied The 

dependence of the two parameters on the char- 
acter of <t>{T) is significant. The mean rate of 
disturbances computed in Paper 1 from equations 
of the impact hypothesis, certain physical assump- 
tions and dynamical evidence comes out 0.25 per 
100 AU (solid line in fig. 3 of paper 1), which is 
in ordcr-of-magnitude agreement with the data 
of table 4. 


PHASE EFFECT AND ALBEDO OF 
COMETARY NUCLEI 

The dimensions and mass of an average ^erratic^ 
comet are important for the impact hypothesis 
for two reasons : 

(1) The impact rate is proportional to the 
collisional cross-section of the comet's nucleus. 

(2) If the impulsive velocity is known, the 
mass of the nucleus determines the magnitude of 
the impulse exerted by a boulder impact, which 
equals the momentum gained by the material 
expelled from the nucleus. The momentum, in 
turn, determines the mass of the boulder. 

An upper limit for a cometary radius can be 
derived from dynamical considerations of a 
cometary splitting. For P/Biela (classed as an 
‘erratic' comet) the requirement that the separa- 
tion velocity be higher than the velocity of escape 
from the surface of the nucleus of radius R and 
mean density p gives a condition 

J?p^/2<1.3 km(g cm~®)^^® (35) 

For p~0.2 g cm”® we find <3 km. If we require 
that the separation velocity exceeds the escape 
velocity from the sphere of action of the nucleus 
we must use another formula (Sekanina, 1968) 
and get for P/Biela 

12 km(g cm”®)^^® (36) 

or 7? <20 km for a low-density snowball. These 
estimates are too crude to be used for the cal- 
culation of an ‘erratic' comet's mass. 

The photometry of faint cometary images at 
large solar distances appears to be more fruitful 


Table 4. 'Erratic' Comets: Total Number and Characteristic Disturbance Rate as a Function of the Disturbance Frequency 

Distribution 


Distribution <^(r) assumed 

Characteristic disturbance 
rate Tc 
(per 100 AU) 

Total of expected 
‘erratic’ comets Nc 

Data fit 

Three dimensional Maxwellian 

0.77 

8 

poor 

Two dimensional Maxwellian 

0.70 

10 

good for A’e^>2 

Damped exponential in T 

0.50 

15 

good for Nc^>l 
good 

Damped exponential in \/t 

0.06 

50 
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for analyzing the dimensions of cometary nuclei. 
However, the practical solution of the problem is 
very delicate not only because of the obvious 
observational difficulties (very faint images; coma 
contamination must be reduced as much as 
possible), but also because of an ambiguity in 
interpretation. 

Disregarding the sources of periodic or quasi- 
periodic variations in brightness (such as the 
shape of the nucleus), we shall deal with the 
geometrical albedo and phase law, the two 
quantities that are determined by the optical 
properties of the nuclear surface and entet the 
reduction photometric formula, from which the 
nuclear diameter is computed. 

There is no chance to obtain direct information 
about the nuclear reflectivity from ground-based 
observations. However, it might be possible to 
determine the phase law from very accurate 
photometric observations. 

The two candidates for the surface texture to be 
considered in reference to the impact hypothesis 
proposed in Paper 1, namely snow of H 2 O and 
asteroid-like compact but porous material, differ 
considerably from each other in both the reflec- 
tivity and the phase variations (Sekanina, 1971), 
A smooth surface of unpacked H 2 O snow has a 
geometrical albedo 0.5, a phase coefficient 
mag deg“^ for small phase angles and generally 
resembles a Lambert surface (Veverka, 1970), On 
the other hand, a typical geometrical albedo for 
asteroids is about 0.15, and the phase coefficient is 
characteristically /3^0.03 mag deg~^ For Icarus, 
this law is still correct at phase angles as large as 
100° (Gehrels et al., 1970). 

Incorporation of the significant phase effect into 
the photometric formula brings the absolute 
brightness up by 0.5"" at phase angle 18°, by at 
42°, and by more than 1"" ever 3 ^ where between 42° 
and 127°, as compared to the Lambert law. 
Because a strong phase effect also implies a lower 
albedo, hence a larger cross section, an average 
asteroid-like nucleus would be larger in diameter 
than a snow covered nucleus by a factor 2.5, 3 
and 3.5, while both nuclei have equal apparent 
magnitudes under equal geometrical conditions at 
phase angles 25°, 48° and 75°, respectively. 

Obviousl}^, the discrimination of cometary 
nuclei by the phase effect can significantly^ improve 
the accuracy^ of the photometric determinations of 


cometary radii, and thus bring down the un- 
certainty in the mass of individual comets by at 
least one order of magnitude. 

Unfortunately, nuclear magnitudes of the 
quality required by the suggested phase-dis- 
crimination method are not available. To illustrate 
the difficulties encountered in an attempt to 
detect phase variations in published sets of 
magnitudes we have compiled table 5 from the 
homogeneous series of photographic magnitudes 
of comets at large solar distances, obtained by 
Roemer and her collaborators (Roemer, 1965, 
1967, 1968; Roemer and Lloyd, 1966; Roemer 
et ah, 1966). The table lists a sample of more 
extensively observed comets of the two types 
considered in Paper 1 (^erratic’ and core-mantle). 
The type identification based on the dymamical 
evidence is given in column 2, the degree of con- 
sistency with the phase-law evidence is com- 
mented on in the last column. Two phase laws 
have been tested in terms of the dispersion in the 
absolute magnitude, Hq: j 3 = 0.00 mag deg“^, an 
approximation good for a smooth snow surface 
and therefore presumably suitable for the 'erratic' 
comets; and ^3 = 0.03 mag dcg~^, which is assumed 
to work reasonably well for the core comets. The 
absolute magnitudes have been computed from 
apparent magnitudes, applying the inverse square 
reduction law. Only magnitudes from solar dis- 
tances larger than 1.2 AU have been made use of, 
so that phase angles have been conveniently^ kept 
within 50°. Table 5 reveals that except for 
P/Arend-Rigaux and perhaps P/Encke, the dis- 
persion in Hq is rather high and the results there- 
fore inconclusive. A particularly bothering trouble 
is a sy^stematic difference between the absolute 
magnitudes at successive apparitions of the same 
comet. This effect is most noticeable in the case 
of P/ Giacobini-Zinner. 

For the sake of comparison we have also cal- 
culated the quantities of columns 3 to 5 of table 5 
for two minor planets of the Apollo type. For 
Adonis, using photographic magnitudes by six 
observers, we obtain respectively +0.7™, dzO.47™ 
and d=0.39™ (the minimum being ±0.36™ for 
^ = 0.0b mag deg-^. For 1960 UA, Object Giclas, 
^\o have used photographic magnitudes by Roemer 
(1965) and photoelectric B magnitudes (reduced 
to the photographic system) by Rakos (1960), 
and found, respectively, +0.5™, ±0.37™ and 
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Table 5.— Phase effect Phase Coefficient {mag /deg)] 


Comet 

Interpretation 
of dynamical 
evidence* 

Average effect of phase 
in absolute magnitude: 
= 0.00)- 
ffoO = 0.03) 
(mag) 

Dispersion in Ho 

Agreement : 
phase effect 
vs dynamics 

P = 0.00 
(mag) 

p = 0.03 
(mag) 

P/Encke 

P / Giacobini-Zinner 
P/Schaumasse 
P/Tempel 2 
P/Arend-Rigaux 
P/Forbes 

P/Schwassmann-Wachmann 2 
P/ Whipple 

core 

erratic 

erratic 

core 

core 

erratic? 

core? 

+ + 
oooooooo 

» if 

OOOOOOOO 

±0.29 
'’0.75 
0.54 ' 
0.60 
0.15 
0.65 
0.59 
0.40 

yes 

yes? 

yes 

no 

yes 

no? 

no? 


“ See: Sekanina (1971), Marsden and Sekanina (1971). ... j ■ 

^ Equals to +0.8™ before perihelion in 1959; and +0.7™ after perihelion in 1959 and in 19oo. 

® Equals to ±0.40™ and ±0.57™, respectively, before perihelion in 1959; and to ±0.34™ and ±0.40™, respectively, 
after perihelion in 1959 and in 1965. 


±0.26"'. This may suggest that a difference of 
O.!*" in the Ho dispersion between the two phase 
laws might already be a meaningful discrimination 
level, if the dispersion itself is within, say, ±0.5"^. 
To obtain more convincing results precise pho- 
tometry must be applied. 

FINAL REMARKS 

We conclude that the hj^jothesis of fairly 
frequent collisions of interplanetary boulders with 
cometary nuclei, suggested in Paper 1 and 
examined from specific viewpoints in the present 
paper, looks reasonably consistent with the 
limited information available on the character 
and rate of disturbances observed in the motions 
of the 'erratic' short-period comets. Precise 
photometry of cometary nuclei, if conducted in 
the future, is believed to improve significantly our 
knowledge of the amount of mass and (mergy 
involved in the sort of collisions under con- 
sideration. We do not exclude a possibility of a 


different interpretation of the observed phe- 
nomena, but we do not see any at present that 
could compete with the impact hypothesis. We 
also feel that collisional processes involving fairly 
large objects of the solar system should be sub- 
jected to extensive investigations rather than 
rejected as ad hoc assumptions without seriously 
considering the chances, effects and characteristics 
of the collisional mechanism itself. 
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NOTE 

Since the time of the lAU Colloquium # 13 our under- 
standing of the problems discussed in this paper has 
further advanced. We know of two more “erratic” comets, 
P/Finley and P/Comas SoM (Marsden, B. G., Sekanina, 
Z., and Yeomans, D. K., 1973 Comets and Nongravi- 
tational Forces. V. Astron. J. 78, 211-225). 



EVOLUTIONARY AND PHYSICAL PROPERTIES OF METEOROIDS 


REFERENCES 

Gehrels, T Roemer, E, Taylor, R. C., and Zellner, B. H., 1970. Minor planets and related 
objects. IV . Asteroid (1566) Icarus, Astron. J., 75, 186-195. 

HAR^v^T, M., 1967. The cloud of interplanetary boulders, in The Zodiacal Light and the Inter- 
planetary Medium, edited by J. L. Weinberg, NASA SP-150, Supt. of Documents, US 
Govt. Printing Office, Washington, 307-313. ■ 

Marsden, B. G., 1969. Comets and nongra\ntational forces II, Astron. J., 74, 720-734. 

, 1970, Comets and nohgravitational forces III, Astron. J., 75, 75-84. 

- >^1971. The motion, splitting and present condition of Biela’s comet, Bidl. American Astron 
Soc., 3, 269. 

Marsden, B. G., and Sekanina, Z., 1971. Comets and npngravitational forces IV, Astron. J. 
76, 1135-1151. * 

Rakos, K. D., 1960. Light variations of the fast moving minor planet discovered on October 22 
1960, by H. L. Giolas, Lowell Obs. Bull, 5, 28-29. ’ 

Roemer, E., 1965. Observations of comets and minor planets, Astron. J., 70, 397-402. 

, 1967. Comet notes, Publ. Astron. Soe. Pacific, 79, 272-273, 372-374, 503-505, 632-634. 

— , 1968. Comet notes, Publ Astron. Soc. Pacific, 80, 106-108, M5-347. 

Roemer, E., and Lloyd, R. E., 1966. Observations of comets, minor planets, and satellites, 

Astron, 71, 443-457. • 

Roemer, E., Thomas, M., and Lloyd, R. E., 1966. Observations of comets, minor planets, and 

Jupiter VIII, Astron. 71, 592-601. ' 

Sekanina, Z,, 1968. Motion, splitting and photometiy of Comet Wirtanen 1957 VI, BuH. Astron 
/nst. Csecli., 19, 153-170. 

, 1971. A core-mantle model for comctary nuclei and asteroids of possible cometary origin, 
in Physical Studies of Minor Planets, edited by T. Gehrels, NASA SP-267, Supt. of Docu- 
ments, U.S. Govt. Printing Office, Washington, 423-428 

Veverka, j. F., 1970. Photometric and polarimetric studies of minor planets and satellites, PhD. 
dissertation, Harvard Univ., 97-107. 

Yeomans, D. K., 1972. Nongravitational forces in comet Giacobini-Zinner, in The Motion, Evalu- 
alion of Orbits and Origins of Comets, edited by G. A. Chebotarev, E. I. Kazimirchak- 
Polonskaya and B. G. Marsden, D. Reidel Publ. Co., Dordrecht, Holland, 181-186. 



